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We investigate the short distance behavior of the nucleon-nucleon (NN) potential defined through 
the Bethe-Salpeter wave function, by perturbatively calculating anomalous dimensions of 6-quark 
operators in QCD. Thanks to the asymptotic freedom of QCD, the 1-loop estimations give exact 
results for the potential in the zero distance limit. We show that the chiral symmetry of the 
gauge interaction implies the existence of an operator whose anomalous dimension is zero for 
a given quantum number. Furthermore we find that non-zero anomalous dimensions of other 
operators are all negative. These results predict the functional form of the NN potential at short 
distance, which is a little weaker than r~ 2 . On the other hand, the computation of the anomalous 
dimension spectrum alone can not determine whether the potential is repulsive or attractive at 
short distance. An additional analytic non-perturbative analysis suggests that the force at short 
distance is indeed repulsive at low energy as found numerically. Some extensions of the method 
are briefly discussed. 
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1. Introduction 

In a recent paper [[I]], the nucleon-nucleon(NN) potentials have been calculated in lattice QCD 
from the Bethe-Salpeter (BS) wave function through the Schrodinger equation. The results quali- 
tatively resemble phenomenological NN potentials which are widely used in nuclear physics. The 
force at medium to long distance (r > 2 fm) is shown to be attractive. This feature has long well 
been understood in terms of pion and other heavier meson exchanges. At short distance, a char- 
acteristic repulsive core is reproduced by the lattice QCD simulation [jl|]. No simple theoretical 
explanation, however, exists so far for the origin of the repulsive core. 

A hint to understand the repulsive core theoretically has appeared in Ref.[||], where properties 
of the BS wave function of the Ising field theory in 2-dimensions, 

<p(r,0) =/(O|a(x,O)a(O,O)|0,-0} in , r=\x\ (1.1) 

are considered analytically. Here 6 is the rapidity of the one particle state. From the operator 
product expansion (OPE) 

ff(jt,0)a(0,0) ~ G(r)l + cr 3 /V(0) + ---, (1.2) 

r— >0 

the BS wave function at short distance becomes 

<p(r,6) ~ Cr 3/4 sinh(0) + <9(r 7/4 ), (1.3) 

which predicts the short distance behavior of the potential as 

(p"(r,G) + sinh 2 G(p(r,G) 3 1 
(p(r,6) ~ ~16V 

The OPE in this case predicts not only the r~ 2 behavior of the potential at short distance but also 
its coefficient —3/16. The potential at short distance does not depend on the energy (rapidity) of 
the state; it is universal. 

In this report, the OPE analysis is applied to QCD, with the aim to theoretically better under- 
stand the repulsive core of the NN potential. 



Mr) = ^777n ~ -TZ3- ( L4 ) 



2. Operator Product Expansion and potentials at short distance 

2.1 General argument 

Let us generalize the argument in the introduction, which relates the operator product expan- 
sion(OPE) to the short distance behavior of the potential. We write the equal time Bethe-Salpeter 
(BS) wave function as 

<PL(J) = (0\O A (r,0)O B (0,0)\E} (2.1) 

where \E) is the eigen-state of the system with the energy E, and Oa,Ob are some operators of the 
system. Here we suppress further quantum numbers of the state \E) other than E for simplicity. 
The OPE of Oa and Ob is written as 

<9 A (F,0)<9 B (0,0) ~ £D^(r)O c (0,0), r= \?\ -0. (2.2) 

c 
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We here assume that the coefficient function D C AB behaves as 

D c AB {r) ~ (-logr)^ (2.3) 

if Oc has the same mass dimension of OaOb- Therefore the BS wave function becomes 

* £(-logr)fe c (£), (2.4) 
c 

whereD c (£) = (0|O c (0,0)|£). 

In the approach of Ref.ljTJ, Bp, the potential is defined by the BS wave function through the 
Schrodinger equation as 

where /X is the reduced mass. Since 

( 2 . 6) 

for the state with angular momentum L, the dominant contribution of the potential for non-zero L 
at short distance is trivially given by 

,«,)*-*£+!> + .... (2 . 7) 

As a non-trivial case let us consider the L = case. There are two cases, depending on the maxi- 
mum value of jG^g defined as j8x = max fi AB . 

(1) j3x 7^ 0: In this case, the potential at short distance is universally given by 

V{r)~- (2.8) 
r l {— logr) 

If j6x > the interaction is attractive at short distance, while it is repulsive if fix < 0. 

(2) px = 0: We denote /3y < is the second largest among j3f B 's. The potential at short distance 
becomes 

V(r)~£lW (29 ) 

which is attractive for Dy(E)/Dx(E) < while repulsive for Dy(E)/Dx{E) > 0. 
2.2 OPE in QCD 

Since QCD is an asymptotically free theory, the 1-loop calculation for anomalous dimensions 
becomes exact at short distance. The OPE in QCD is written as 

<9 A (r,0)<9 B (0,0) = Y, D AB(r,g, m, \X )O c (0,0) (2.10) 
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where g (m) is the renormalized coupling constant (quark mass) at scale jj,. In the limit that r = 
e *R — > ( t — > oo with fixed R), the renormalization group analysis leads to 

D c AB (r,g,m,Ll) ^ (-l^gHogr/^/^MBiRAO^), (2.11) 

r— >0 

where 

'■"-tsfO'- 2 ?) 

is the QCD 1-loop beta-function coefficient, and 

C,(l) (1) , (1) (1) ^ io\ 

Here is the 1-loop anomalous dimension of the operator Ox- An appearance of D^ B (/?,0,0,^i) 
on the right-hand side tells us that it is enough to know the OPE only at tree level. From the above 
expression, j6f B in the previous subsection is given by 

y c,(i) 

= pry- ^ 

Therefore our task is to calculate y^ 1 ' for 3 and 6-quark operators. 

3. Anomalous dimensions for 6-quark operators 

3.1 6 quark operators 

The OPE of two baryon operators at tree level is given by 

B 1 {x)B 2 {0) = B l (0)B 2 (0)+i"(^ 1 (0))B 2 (0) + V/(^ v fi 1 (0))8 2 (0) + - (3.1) 

where the first term corresponds to the L = contribution, the second to the L = 1 , and so on. In 
this report we consider the L = case (the first term) only. We denote the general form of a gauge 
invariant 3-quark operator as 

= B rM = e ahc q a a f (x)q h f{x)q c 7 \x) (3.2) 

where a,j3,y are spinor, f,g,h are flavor, a,b,c are color indices of quark field q. The 6-quark 
operator is constructed from two 3-quark operators as 

B F T ^ 2 {x)=B F T \(x)B%(x) (3.3) 

where = a,/?,); and F; = figihi (i = 1,2). Since quarks are fermions, there are some linear depen- 
dencies among 6-quark operators. We have to determine a set of independent 6-quark operators. 
It is not so easy, however, to find them using a quark field basis. Instead we use a simpler method 
mentioned below. 

As the choice of the gauge fixing in perturbation theory, we take the covariant gauge with 
gauge parameter A. Since both the 3-quark operator B^ and 6-quark operator S^'^? are gauge 
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invariant, Jab must be independent of the gauge parameter X. The A dependent term in the 

c (11 

calculation of y A £ at 1-loop becomes 

+ * t^$&)> (3-4) 

where the i-th index of abc and the j-th index of def is interchanged in (abc, def) [ij] . For example, 
(ri,r 2 )[ll] = «2j8i7i,aij8 2 72 or (ri,r 2 )[21] = aia 2 7i,j8ij3 2 72. Note that the interchange occurs 
simultaneously for both ri,r 2 and F\,F% in the above formula. The gauge invariance implies 
eq.(fO|) = 0, which yields constraints among the 6-quark operators. For example, let us consider 
the case that ri,T 2 = aa/3, a/3/3 and Fi,F 2 = ffg,ffg, for which the constraint becomes 

■3 I \offgJfg I fl _r)\T>ffgjfg J_T>fffJgg I /O _ l\l>fgg,fff \ _ A 

3 ^ace/3,cej3j3 ^ \ 3 z ^aa/3,ajSj3 ~ rD aaa,ppp V z l ! D app,aap) ~ u 

^ADffgJfg J_Dfff,fgg j_T>fggJff _A , T Ci 

^ ^ D aap,aPP ^ D aaa,ppp ~ r °app,aap ~ u > ^"^ 

where minus signs in the first line come from the property that B^yr = — Bp*' p . There are no 
further relations among 6-quark operators beyond (JT4|). 

3.2 1-loop contributions 

Only the divergent part of the gauge invariant contribution at 1-loop is necessary to calculate 
the anomalous dimension of 6-quark operators at 1-loop. The building block of 1-loop calculations 
is the gluon exchange between two quark lines. If both two quark lines belong to one operator, 
either fip or B[? , the contribution is canceled by the renormalization factor of the 3-quark operator 

Bp, so that the divergent term does not contribute to 7^g = Ya^ + Yb^ ~ Yc^ • ^ one 1 uar k l me 
comes from Bp' ( and the other from B^?, the divergent term contributes to 7^ . Suppose that one 
quark line has indices ((Xaj/a) at one end and (ofi,/i) at the other end and the other quark line has 
(a B ,/ B ) and (a 2 ,/ 2 ). The divergent contribution from the 1-gluon exchange can be expressed as 

g 2 1 

-^~2 1 { 5 fifA 5 fzfB [Sa iaA 8a 2 a B ~ 28 a2 a A 8 ai a B ] + 38f 2 f A S fl f B [8 a2 a A 8 ai a B ~ 2S ai a A 8a 2 a B ] } (3.6) 

for (ai , a 2 ) € (R,R) or (ai , a 2 ) € (L,L), where R and L means the right and the left handed compo- 
nent of the spinor indices. Other combinations, (R,L) or (L,R), vanish. This property comes from 
the fact that the gluon coupling to quarks is chirally symmetric and the chirally non-symmetric 
quark mass term does not contribute to the divergence. 

3.3 Chiral decomposition of 6-quark operators 

The physical nucleon operators are constructed from general 3-quark operators as 

B f a=B^ y (Cy 5 ) fiy (iz 1 y h (3.7) 

where C is the charge conjugation matrix, f,g,h are u or d, and T 2 is the Pauli matrix in the flavor 
space. The spinor index a is restricted to the positive energy component such that a = 1 , 2 in the 
Dirac representation of the y matrices. 
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In this report we consider L = two nucleon states, which are l So and 3 Si. Here we use the 
notation 2S+l Lj where S is the total spin, L is the orbital angular momentum and J is the total 
angular momentum. The 6-quark operator for 1 So, which is the spin-singlet and isospin-triplet 
state, and for 3 S\ (the spin-triplet and isospin-singlet state) are given by 

BB(%) = (ia 2 ) ap B f a B f p , flfl( 3 5i) = {ix 2 ) fg B f a B 8 a , (3.8) 

where the summation is taken for the repeated index. Both 6-quark operators have the following 
chiral decomposition: 



BB — BliBh + BiiBir + BuBrl + BlrBlr + BuBrr + (L^R) 
where Bxy means B a wy\ with a € X and [/$, y] G Y for X,Y = R or L. 
3.4 Anomalous dimensions 



(3-9) 



We now give our main results in this report. We define 



c,(i) 
Tab 



(i) , (i) (i) 
Ta +Tb -Tc 



32k 2 



y- 



(3.10) 



The eigen-operators of the anomalous dimension matrix y are found to correspond to the chirally 
decomposed operators in the previous subsection. We give the eigenvalue of each operator in 
table [lj which shows that the operator with zero anomalous dimension always exists and other 
anomalous dimensions are all negative for both and 3 S\ states. This corresponds to the case (2) 
of the general discussion in the section \2.l\ 

(2) (f3 x = 0,p Y = -3/(33 - 2N f )) for % and (/fe = 0, j8 F = -1/(33 - 2N f )) for 3 Sl 

The appearance of zero eigenvalues in both 1 Sq and 3 S\ states can be understood as follows. 



As mentioned in section |3.2| , 1-loop contributions to the y A ^ exist only if spinor indices of two 
quark lines, one from T\ the other from 1^2 in By^By 2 , belong to the same chirality (left or right). 
Since BuBrr + BrrBh has no such combination, y A ^ is always zero for this type of operators. 
As pointed out before, this property is the consequence of the chiral symmetry in QCD interactions. 



Table 1: The eigenvalue y for anomalous dimension of each eigen operator in 1 So and 3 S\ states. In the table 
(XY,ZW) means B XY B Z w + {R^L). 



7(%) 
7( 3 Si) 



(LL,LL) (LL,LR) (LL,RL) (LR,LR) (LR,RL) (LL,RR) 



-12 

-28/3 



-4 
-4/3 



-16/3 








4. Conclusion 

The OPE and renormalization group analysis in QCD predicts the universal functional form 
of the nucleon-nucleon potential at short distance: 

V f r ) ~ Dy ( £ ) Zh r ^o (4i) 

[ > ~ D x (£)r2(-logr)i-fr' U ' 
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which is a little weaker than a 1 /r singularity. We obtain 

fr^—srW Ww—srW <42) 

The anomalous dimension spectrum, however, cannot alone tell whether the potential at short 
distance is repulsive or attractive. If we evaluate Dx(E) and Dy{E) by the chiral effective theory 
at the leading order {i.e. the tree level), we obtain 

DAE), 2E Dy(E) /3c 2m N 

r , / j-r \ { ?>o) = , t-7m =ir, ( 4 -3) 

D X (E) m N D X (E) E 



where E = y p 1 + nv\ and p is the relative momentum of two nucleons. We find a repulsive core 
for both states. In particular, in the low energy region such that p 2 <C m^, the repulsive potential at 
short distance is almost energy independent, since 

Furthermore no low energy constant of the effective theory is required at leading order to obtain 
the above result. 

There are several interesting extensions of the analysis using the OPE. The extension of the 
OPE analysis to the 3-flavor case may reveal the nature of the repulsive core in the baryon-baryon 
potentials. Since quark mass can be neglected in this OPE analysis, the calculation can be done in 
the exact SU(3) symmetric limit. It is also interesting to investigate the existence or the absence of 
the repulsive core in the 3-body nucleon potential. In this case, we have to calculate anomalous 
dimensions of 9-quark operators at 2-loop level. More precise evaluations of the matrix element 
Dx{E) = (0\Ox\E) would be preferable. The most straightforward extension is to analyze the 
tensor force and LS force by the OPE. Preliminary results indicate that 

V T (r) ~ Co + di-logr)^- 1 , V LS (r) ~ (4.5) 

where LS force has the strong attractive core at short distance, while no repulsive core exists for 
the tensor potential. Absence of the repulsive core predicted in the tensor force is consistent with 
recent numerical simulations g]. 
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